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ABSTRACT 

Runaway stars produce shocks when passing through interstellar medium at 
supersonic velocities. Bow shocks have been detected in the mid-infrared for 
several high-mass runaway stars and in radio waves for one star. Theoretical 
models predict the production of high-energy photons by non-thermal radiative 
processes in a number sufficiently large to be detected in X-rays. To date, no 
stellar bow shock has been detected at such energies. We present the first de- 
tection of X-ray emission from a bow shock produced by a runaway star. The 
star is AE Aur, which was likely expelled from its birth place by the encounter of 
two massive binary systems and now is passing through the dense nebula IC 405. 
The X-ray emission from the bow shock is detected at 30" to the northeast of 
the star, coinciding with an enhancement in the density of the nebula. From the 
analysis of the observed X-ray spectrum of the source and our theoretical emis- 
sion model, we confirm that the X-ray emission is produced mainly by inverse 
Compton up-scattering of infrared photons from dust in the shock front. 
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Introduction 



Runaway stars move throu gh space with an extremely high velocity relative to the 
surrounding interstellar medium ( lBlaauwlll96ll ). A proposed scenario for the ejection of stars 
at high velocities from their birth place is close interaction between two binary systems (e.g. 
Hoogerwerf et al.l 120001 ) . This process may produce ejection velocities close to the escape 
velocity at the surface of the most massive star in the system: up to 100 — 1000 kms -1 
( iLeonardl Il99ll ; iGvaramadze fc Bomansl 120081 ) . The enc ounter of the proper mo tions of two 
or more stars favors the dynamical ejection scenario ( Hoogerwerf et al.l 120001 ) . A typical 
example of a system of runaway stars escaping from a close interaction of two high -mass 
binaries is the group formed by AE Aur, \x Col and i Ori (e.g. iGualandris et al.l 12004 ). 



The passag* 


3 of a 


Mohamed et al. 
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[e.g. 



were i dentified in the past with the Infrared Astronomica l Satellite (IRAS; Ivan Buren fc McCray 



1988; van Buren et al. 



1995 



Noriega-Crespo et al.lll9970 and, more rece ntly, with the Spitzer 



space telescope ( Kobulnicky et al. 201ol ) and other infrared missions (e.g. Comeron fc Pasquali 
2007h . 



Recently, iBenaglia et al.l (120101 ) detected non-thermal radio emission in the massive 
runaway star BD+43° 3654. This emission was associated with cooling of energetic electrons 
through synchrotron radiation and the authors hypothesized that some bow shocks may 
also produce h igh-energy emission that m ight be detectable. This problem was treated 
theoretically in Idel Valle fc Romero! ( 120121 ). who predicted detectable X-ray and gamma-ray 
emission associated with particle acceleration at the b ow shock. No o bsery ational constraints 
were found to this phenomenon to date. Very recently, iTerada et al. J2012[ ) reported the non- 
detection of X-ray emission from the bow shock produced by BD+43° 3654 in a long time 
exposure X-ray observation with Suzaku. In this letter, we report the first detection of 
hard-energy emission from a bow shock. 



2. The Runaway Star AE Aur (HIP 24575) 



AE Aur was ejected from its birth place in the Orion nebula cluster around three 
million years ago, very probably by the encounter of two massive binary systems. As a 
result, the stars AE Aur and \i Col (both with spectral type 09.5) were expelled at high 
yelocities, while the syst em i Ori remained as a bounded binary system with high eccentricity 
( IHoogerwerf et al.ll2000l ). 



AE Aur has recently encountered IC 405, a dense molecular nebula with Hydrogen den- 
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sity riYL ~ 3 cm -3 (IPeri et al.ll2012l ). As a result of this enc ounter, a bow shock was prod uced. 
The bow shock was detected in t h e mid -infrared band by Ivan Buren fc McCrayl (119881 ) with 
IRAS and later by France et al.l (120071 ) using Spitzer. Some important para meters of the 



star for our study are the terminal velocity o f the wind 
201lh . the mass-loss ra te (M ~ 10~ 7 M ^yr" 1 ; 



1500 kms , iHubrig et al. 



Fullerton et al.ll200q ). and the velocity of the 



star (u* f» 150 kms 1 ; IPeri et al.l 120121 ) . In Tabled! we list a set of physical parameters of 



AE Aur, including some assumptions on the characteristics of relativistic particles adopted 
for our model. 



3. X-ray Observation and Data Analysis 

For our work, we analyzed the archive XMM-Newton EPIC observation of AE Aur (ID 
0206360101, PI F. Damiani), performed on 2004 September 10 in Full Frame Mode with 
the Medium filter. The data were processed with the Science Analysis System (SAS V12.0). 
Light curves, images, and spectra, were created by selecting events with PATTERN<12 for 
the MOS cameras, PATTERN<4 for the pn camera, and FLAG=0 for both. To reduce 
the contamination by soft proton flares, the original event files were screened by using the 
sigma-clipping algorithm (ESPFILT tasks). The screened MOSl/MOS2/pn exposure times 
are 52/54/38 ks, respectively. Spectral analysis was performed in the 0.3 — 7 keV energy band 
using XPEC VI 2. 7 (the Ancillary Response Files were generated with the SAS ARFGEN 
task). MOS and pn spectra were rebinned to achieve a signal-to-noise ratio per bin > 5a. 
All the reported errors are at the 90% confidence level. 



Figure [T] shows the WISE 12.1/xm image (in red, see also IPeri et al.l l2012h of the bow 
shock of AE Aur together with the corresponding EPIC pn count image in the 1 — 8 keV 
band (left panel, in green) and the pn median photon energy map, i.e., an image where each 
pixel holds the median energy of the detected pn photons in the 0.3 — 8 keV band (right 



pane l, in green; for details on the procedure adopted to produce the map, see iMiceli et al. 



20081 ). A bright X-ray source, ~ 30" northwest of the star is visible in the left panel of Fig. 
[T] (BS region) and appears embedded in the infrared bow-shock. We verified that BS has 
no point-like counterparts in the optical and infrared bands. The median photon energy 
map clearly shows that the X-ray emission from the BS source is significantly harder than 
that of AE AmQ. Namely, the average value of median photon energy is ~ 850 eV in region 
BS and ~ 750 eV in region AE Aur. Moreover, this map reveals a bow shaped structure 



1 We verified that the MOS cameras provide consistent results. Here we present pn images because of the 
better statistics 
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Fig. 1. — Left panel: WISE 12.1//m image (red) and EPIC pn count image of AE Aur in the 
1 — 8 keV band (green). The bin size of the X-ray image is 4". The regions selected for the 
spectral analysis and the proper motion vector are superimposed. BKG is the annular region 
from which the background spectrum was subtracted (see text). North is up and East is to 
the left. Right panel: same as left panel with pn median photon energy map in the 0.3 — 8 
keV in green. In the median photon energy map each pixel holds the median energy of the 
detected pn photons band. Only pixels with more than 4 counts we re taken into accou nt 
and the image is smoothed with a = 8", according to the formulae in lMiceli et ali ( 2008 ). 



(not directly visible in the count image) characterized by hard X-ray emission and somehow 
reminiscent of the infrared bow shock. 




Energy [keV] 



Fig. 2. — Background-subtracted EPIC pn spectrum of region BS of Fig. [JJ with the theo- 
retical model described in Sect JH and the corresponding residuals (Observed - Model). 

We first analyzed the stellar spectrum by extracting the pn and MOS spectra from 
region "AE Aur" of Fig [TJ We subtracted to the source spectrum a background spectrum 
extracted from a nearby region (out of the field of view of Fig[T]). The fittings were performed 
simultaneously on both MOS spectra and on the pn spectrum by adopting an absorbed APEC 



model (optically thin coronal plasma in collisional ionization equilibrium, ISmith et al.ll200l[ ) 
based on the 2.0 release of the AtomDB database. The interstellar absorption was described 
by the PHABS model in XSPEC. The stellar spectrum can be fitted (reduced x 2 = 1-32, 



with 315 d. o. f.) by an isothermal plasma with temperature kT = 0.225_ 004 keV, emission 

0.5 

o.i 



measure EM = 2.8+04 x 10 55 cm 3 (assuming a distance D = 550 pc, in agreement with 



Peri et al. 2012 ). and oxygen abundanc^l 0/0 Q = 0.26 ± 0.02. The best-fit value of the 
column density is n# = 4.2 ± 0.2 x 10 21 cm -2 and the unabsorbed X-ray luminosity in the 
0.3 - 10 keV band is L x = 1.9 x 10 32 ergs" 1 . 

We then analyzed the X-ray spectrum of the BS source (the extraction region is shown 
in Fig. [1]). The contamination from AE Aur photons scattered by the telescope PSF in the 
BS region was removed by extracting the background spectrum from the region BKG of Fig. 
[TJ This region consists of an annulus centered at the star and with inner and outer radii that 



2 We verified that the plama metallicity is ~ 0.35 and we fixed the abundances of all other elements to 
this value. 
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were chosen to contain pixels at the same distance from the star than the bow shoclo Figure 
[2] shows the pn spectrum of region BS. Only 180 pn counts (943 pharcmin -2 ) were detected 
in this region, while the spectra from the MOS cameras do not provide useful constraints 
and are not presented here. The BS spectrum can be well fitted by an absorbed power law 
model with photon index T = 2.6i 5 an d normalization iV = 6 ± 2 x 10~ 6 photons keV -1 
cm -2 s _1 at 1 keV. Since the absorbing column is consistent with that found for AE Aur, we 
fixed it to % = 4.2 ±0.2 x 10 21 cm -2 in the fittings, to obtain tighter constraints on the other 
model parameters. We used Cash instead of \ 2 Statistics for the fitting, that is recommended 
when the number of counts in a channel is small. We used the task GOODNESS with 1000 
Montecarlo simulations to determine the goodness of our fit. The value obtained for our fit 
is 2.3% of realizations with less best fit statistic than that obtained with the fitted model 
(i.e. 97.7% confidence interval). 

The BS spectrum can also be well fitted by a very hot thermal component with kT = 
2.41q 8 keV and EM = 4 ± 1 x 10 52 cm~3 (assuming that BS is at 550 pc). Such extremely 
high temperature cannot be associated with the stellar wind, as shown below. If we assume 
that the hot plasma fills the whole region, we derive from the EM a lower limit for the 
particle density, n ~ 3 cm -3 , that is much higher than that expected for a stellar wincE). 
Temperatures higher than 1.6 keV can be indicative of a cataclysmic variable star or of a 
very energetic coronal flare, though we verified that the BS light curve is consistent with a 
constant source. Given the low count-rate of the source, statistical fluctuations are quite large 
and deeper exposure are necessary to obtain a proper study of the light curve. In conclusion, 
present data do not allow us to discriminate between a thermal and a non-thermal scenario 
for the BS X-ray emission and deeper observations are necessary. Nevertheless, the lack 
of infrared and optical counterparts and the extremely high temperature derived from the 
thermal model make the association of BS with a foreground/background stellar object less 
likely. On the contrary, the spatial correlation between BS and the infrared bow shock, 
the bow-shaped morphology of the hard X-ray emission around AE Aur, and the spectral 
properties of BS (see next section) concur in suggesting that BS is indeed the result of 
non-thermal emission originating at the bow shock of the runaway star AE Aur. 



3 Pixels inside the BS region were extracted from the background extraction region (BKG). 
4 If the plasma does not fill the entire region we get even higher density values. 
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4. Non-Thermal Radiative Model 

In order to explain the observed X-ray emission from the bo w shock of AE Aur, we have 



applied the non-thermal radiative model of £ Oph, developed by Idel Valle fc Romero! (120121 ) 



with the parameters given in Table [T] (see als o § El) . By adop ting stellar parameters of an 



09.5V star and a ISM medium Nr = 3 cm 3 ( IPeri et al.ll2012[ ). we obtain that the runaway 



star produces a bow shock in the interstellar medium with a stagnation point located at 
approximately 17.1 x 10 3 AU. This estimation agree very well with the value measured from 
the X-ray observation (~ 17 x 10 3 AU). The reverse shock in the stellar wind is adiabatic 
and provides an efficient site for particle acceleration through a Fermi I type acceleration 
process. Since the gas is compressible, we assume that it is in su b-equipartition, getting an 



estimate of ~ 10 4 G for the magnetic field in this medium (see Idel Valle fc Romero! 12012 
for details). 

The available power in the stellar wind is ~ 4 x 10 33 ergs -1 . Assuming that only a small 
fraction of ~ 1 % is transformed in relativistic particles (we assume equal power in leptons 
and hadrons), we have been able to reproduce the observed X-ray spectrum with a power- 
law electron injection of index —2.6, through inverse Comptom up-scattering of IR photons 
coming from dust that is heated by the forward shock and the thermal radiation. This model 
produces a good fit to the observed spectrum of the bow shock (see Fig. [2]), with reduced 
X 2 = 0.43 with 5 d.o.f. (see § [2]). Models with higher ratio of proton to electron power are 
possible. Since the cooling time of the protons is much longer, they are removed from the 
radiation zone without having time to make a significant contribution to the spectral energy 
distribution. So, in order to keep the goodness of the fit, a larger total power is required for 
proton-dominated models. In a case similar to that of the Galactic cosmic rays, where the 
mentioned ratio is 100 instead of 1, we should consider that about 0.3 of the available power 
goes to relativistic particles. This requires very efficient acceleration. The case with a = 1 is 
energetically more conservative, but we cannot rule out models with a up to 100. In all cases, 
the shape of the SED is basically the same and determined by the electron population. 

In Figure El we present the broadband spectral energy distribution implied by our model. 
Inverse Compton losses dominate, producing a steepening of the non-thermal spectrum that 
results in the soft power-law X-ray emission observed. The dominance of the radiative 
losses requires a low level of convection of the relativistic particles that can result from the 
development of Kelvin-Helm holtz instabilities in the bow shock (see conditions for these 



instabilities in bow shocks in lAraudo et al.ll2010l ). 



The SED predicted for the runaway star has little and steep synchrotron radio emission 
and a peak of around ~ 10 29,5 erg s _1 at soft X-rays. Contrary to the case of ( Oph, 
located closer and with different physical parameters, the expected emission at gamma-rays 



-8 - 




Fig. 3. — Spectral Energy Distribution of non-thermal radiati ve processes for the AE Au r 
bow shock. The model was computed following the work by del Valle fc Romero! ( 2012 ). 
The continuous line represents the inverse Compton process produced by the dust (See § H] 
for details). 



is negligible. Hence, the best test for the proposed model is the radio detection through deep 
interferometric observations at long centimetric wavelengths. 



5. Conclusion 



We detected a hard X-ray source spatially correlated with the infrared bow shock of the 
runaway star AE Aur. The source is consistent with being a point-like steady source with 
non-thermal emission, though deeper observations are necessary to ascertain its morphology, 
its variability (if any), and its photon index. Nevertheless, the analysis of the X-ray source 
at ~ 30" to the northeast of the runaway star AE Aur showed that it is very likely the 
X-ray counterpart of the bow shock produced by the passage of the star through the dense 
molecular nebula IC 405 and that was detected in the infrared with IRAS, Spitzer and WISE. 
We developed a radiative non-thermal model specifically tuned to describe the AE Aur bow 
shock properties and obtained a very good agreement between model and observations. 
According to our model, the X-ray emission would be produced by inverse Compton process 
of accelerated particles in the nebula dust. High-energy emission from bow shocks produced 
by runaway stars has been predicted by theoretical models. This is the first time that X-ray 
emission is detected in one of such bow shocks. Future X-ray and radio observations of this 
source will permit to further constrain its properties. 
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Table 1: Parame ters for the non-thermal radiative model for AE Aur 
Parameter value 



Ro 


Standoff radius 


1.7 xlO 4 AU 


M w 


Wind mass loss rate 


lCT 7 M yr- 1 


a 


Hadron-to-lepton energy ratio 


1 


<?rcl 


Content of relativistic particles 


0.007 


a 


Particle injection index 


2.6 


v w 


Wind velocity 


1.5xl0 8 cm s~ x 


L 


Available power 


4xl0 33 erg s _1 


B 


Magnetic field 


l.lxlO- 4 G 


V, 


Star velocity 


150 km s _1 



